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ABSTRACT 

An experimental study has been made of the feasibility of 
following the precipitation process in a polycrystalline aluminun- 
24% copper alloy. The intermediate structures formed during ageing 
have previously been determined by single crystal methods to pass 
through the following sequence: 

G.p( + G.P[2zi— 6 — 8 

The sequence of structures was followed by analysis of the diffrac-— 
tometer line profile of the &(@00) peaks. The experimental 
profiles were corrected for instrumental broadening by the Stokes 
method. One-dimensional Patterson projections of the electron 
density distribution in the alloys were synthesized from the 
corrected profiles. 

The complete sequence of intermediate precipitates could not 
be realized in this investigation. However it was possible by this 
method to distingish between the G.P. [\] and Q) structures by 


observing the height and shift of the Patterson function origin peaks. 
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ie Introduction. 

The phenomenon of precipitation hardening holds great commer- 
cial interest. It is a method of enhancing the mechanical properties 
of a metal without resorting to other processes such as work hardening, 
dispersion hardening or a martensitic type transformation. It elim- 
inates the concern regarding the hardenability of a metal. Metal 
parts may also be formed or fabricated in a soft condition and then 
hardened by a relatively simple ageing treatment while retaining 
excellent dimensional stability. One example of the importance of 
precipitation hardening is its application to aluminum alloys which 
has resulted in a great number of usable high strength materials. 

The mechanism of precipitation hardening is not at all simple. 
The most basic concept involves the formation of a stable phase 
dispersed throughout a matrix. This does not include dispersion 
hardening, which is due to the presence of a second phase dispersed 
in the matrix. However, it is not a precipitation process. Even 
in the basic concept of precipatation there are constraints imposed 
on the nucleation and growth of ae a phase because it is being 
formed in a solid solution. In general, the process is much more 
complicated than this. It varies with each alloy system but usually 
involves the precipitation of a metastable phase, which later gives 
way to a stable one. There may be as many as three intermediate 
stages before the stable precipitate is formed. The intermediates 
may or may not be true phases, depending on the alloy system. 

The intermediate stages hold great interest. In many alloy 
systems these intermediates have the greatest influence on mechanical 


properties, and the formation of a stable precipitate results in a 


t 





decrease in hardness anc strength. 

The Al-Cu alloy represents an ideal system in which to study 
the age-hardening process. It is representative of many systems for 
it displays both substitutional and displacement disorder during 
ageing. The copper atoms have a greater scattering factor and a 
smaller radius than the aluminum atoms (Feu/Fa, = 29/13 and 2.56 A 
compared to 2.86A ). The difference in atomic size results in the 
displacement disorder while the greater scattering factor of copper 
allows its effects to be detected by x-ray analysis. 

The age-hardening system in Al-Cu has been studied more than 
any other. AS a result its sequence of structures is known in consid- 
erable detail.! 

The mechanism of the decomposition of the quenced supersatur- 
ated solid solution shows four different stages whose relative roles 
depend strongly on the annealing temperature. In general the sequence 
may be represented as follows: 

GP] —~ apizli—~ o— 6 
Two stages, 6 andQ, correspond to true precipitates of composition 
CuAly, where 8 is the equilibrium phase. The G.P. refers to a 
"Guinier-Preston zone", which may be defined as a region rich in 
solute atoms which is formed on certain planes in the solvent matrix 
and retains what is essentially the structure of the matrix. G.P.|\\ 
is a copper rich region of plate-like shape formed on {\Oo} planes 
of the aluminum matrix whose maximum size is a few atomic planes 


2 
thick and up to about 100A in diameter. G.P.(2] displays features 


lJ, M. Sileock, T. J. Heal and H. K. Hardy, J. Inst. Metals 82, 
240-242 (1953). 





of a zone and a true phase. As these platelets reach a thickness 

of about 30 A they become unstable and dissolve. They may be considered 
as having a lattice similar to that of the aluminum matrix. To a first 
approximation it is merely a superlattice of aluminum, tetragonal 

with c/a = 2 where a is the parameter of aluminum. 

The above structures do not necessarily evolve from one another 
in sequence, but two (or more) may occur simultaneously in the same 
specimen. The conditions have been determined in which the three 
metastable phases and the equilibrium phase appear as a function of 
time and temperature for different concentrations of Cu in Al.° 

A variety of x-ray techniques have been used in the numerous 
precipitation investigations. Single crystals have generally been 
employed using both white and monochromatic radiation. Some work 
has been done using normal powder methods. The small angle scatter 
method has been applied with great success by Guinier.? The tech-— 
niques using Single crystals and monochromatic radiation have 
advantages of ease and certainty of interpretation in terms of a 
reciprocal lattice plot, but require sophisticated equipment and 
precise experimental technique. Some of the results obtained by 
the above methods have recently heen confirmed by electron micro- 
scopic photographs. 4 The use of an x-ray diffractometer in the 
analysis of polycrystalline specimens is attractive because of the 


relatively simple procedures involved and the possible extension of 





“Ibid. 9 Ppp. 243-45. 
8a. Guinier, Ann. Phys. 12, 161 (1939). 


4G, Thomas, Transmission Electron Microscopy of Metals, pp. 236-248. 





the technique to precipitation studies in commercial alloys. Hence, 
it is the purpose of this investigation to follow the ageing process 
in a polycrystalline Al-4%Cu alloy by means of a diffraction line 
shape analysis and Patterson projection of the electron density 
distribution within the alloy. 

The data necessary to form a Patterson function is available 
from a diffraction line profile, as is shown below. However, a true 
line profile must first be obtained by correcting the observed profile 
for instrumental broadening, due mainly to the slight spread of wave- 
lengths in the incident radiation. A method has been derived by 
Stokes using the convolution theorem of mathematics to correct the 
experimental line shape without any prior assumptions concerning the 
corrected line shape.° The method consists essentially of obtaining 
the Fourier coefficients for the broadened and reference curves, 
operating on these coefficients (termed unfolding), and then synthe- 
sizing the corrected curve from the unfolded coefficients. This 
technique eliminated broadening due to Ky doubling as well as other 
forms of broadening. 

It has been shown that the true line profile is proportional to 
the Fourier transform of the square of the absolute magnitude of the 
structure factor, \Frcal .© The structure factor Fug, in general a 
complex quantity, represents something real about a crystalline 


material for it is the Fourier transform of the electron density in 


Oa, R. Stokes, Proc. Phys. Soc. London 61, 382 (1948). 


Sh. J.) Wilson, X-Ray Optics, p. 107. 





the material. ! Experimentally only female is determined which yields 
no information concerning its phase. If the phase relations could be 
determined experimentally there would be no problem in finding the 
arrangement of atoms in any crystal, no matter how complicated. Of 

the methods proposed to circumvent the phase problem, the Patterson 
function has emerged as the most powerful and adaptable. This function 


8 If one considers 


is a Fourier series whose coefficients are (eens ‘ 
the Fourier coefficients \Fwteal as derived by forming products of 
each Fourier coefficient Fryer with its complex conjugate, the electron 
density function can be derived directly from the Patterson function: 
The Patterson function is essentially a display of all the vector 
lengths between atoms of a unit cell. All the vectors of one length 
will then appear as a peak in this display. Since the reverse vectors 


are also displnyed, it is necessary to sort out the duplicated peaks 


in order to arrive at the electron density function. 


26 Alloy preparation. 

The Al-4%Cu (by weight) alloy was prepared using 99.99% ingot 
Al and 99.94% Cu shot. The aluminum was first melted in an induction 
heated graphite crucible, and then the copper shot was added. The 
molten metal was held between 700° and 720°C for 45 minutes to ensure 
complete solution of the copper. The alloy was cast into a flat 
rectangular zirconia boat, preheated to 700 C, and allowed to solidify 


at room temperature. The ingot thus formed weighed 125 grams, 


em, PO 


Ta. Guinier, X-Ray Diffraction, p. 92. 


ou J. Buerger, Vector Space, pp. 5-10. 
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measuring approximately 4" x 24" x %". 

Microscopic examination of the ingot revealed a large grained 
structure containing intergranular theta phase and no evidence of 
undissolved copper. 

The ingot was homogenized at 542°C for 24 hours and furnace 
cooled. Grain refinement was accomplished by cold rolling and 
recrystallizing at 600°C for one hour. The amount of reduction 
per pass through the rolls was limited to about 204 due to work 
hardening, so the refinement treatment had to be repeated three 
times, resulting in an ingot .17" thick, and containing equiaxed 
small grains with finely divided theta phase dispersed throughout. 

Eight hardness specimens were cut from the ingot, each being 
about 1.75" x .5" x .17" in size. A hole was drilled in one corner 
of each specimen and top and bottom surfaces were smoothed on no. 0 


and no. 3/0 metallographic paper. 


3. Solution treatment and age hardening. 

The eight specimens were strung on a wire with a $" separation 
to ensure rapid heat transfer both on heating and quenching. They 
were solution treated in a salt bath for six days at 536°C, which 
is about 36° above the solid solubility of 4% copper in aluminum. 2 
Following the solution treatment the specimens were quenched in 
room temperature water, and the as quenched hardness was taken 
immediately on one specimen. The surface irregularities caused by 


the salt bath on the remaining specimens were removed by lightly 


Qn, Hansen, Constitution of Binary Alloys, Pp: 84. 





polishing on no. O and no. 3/0 paper. They were then suspended in 
an oil bath at 165°C for ageing. The as quenched hardness specimen 
was aged at room temperature. 

Microscopic examination of the as quenched hardness specimen 
following solution treatment revealed a small amount of undissolved 
theta phase in the structure. A lattice parameter measurement on 
the &(422) diffraction peak of this specimen indicated that about 
24% of the copper had gone into solution. The small amount of undis- 
solved theta would have no effect other than to reduce the peak 


hardness obtained during ageing. 


4. Diffraction procedure. 

As each specimen reached its respective ageing time, it was 
quenched in room temperature water and hardness readings were taken 
immediately. No less than six readings were taken on both sides to 
ensure a good average. The surfaces were then given a rough electro- 
polish to remove the oxide layer and to ensure that the surface was 
as nearly a true representation of the bulk of the specimen as possible. 
Chemical macro-etching was not used due to its preferential removal of 
copper from the alloy. 

The &(200) peak of each specimen was recorded using a Norelco 
x-ray diffractometer with Ni filtered Cuka radiation. The diffracted 
radiation was counted by a xenon proportional tube in conjunction with 
a pulse height analyzer set to pass the central 90% of the pulse 
distribution. 

The reference peak for the Stokes! unfolding was obtained from 


a recrystallized specimen of the pure ingot aluminum. It was possible 





to get a well resolved Ka doublet for the (200) peak from this specimen. 


De Treatment of data. 

The intensity of each diffraction peak was recorded as a function 
of R@Q. Due to the sharpness of the reference peak it was necessary to 
use an interval of .0125° 28 between points to describe the peak by a 
smooth curve. It was possible to use a point interval of .025° 26 to 
describe the broadened peaks. A total interval of 1.8° was used on 
all broadened peaks except the as quenched and aged one day specimens, 
where it was possible to use slightly shorter intervals due to their 
narrowness. The maximum value of the diffraction peak was centered 
in the total interval in all cases. 

The unfolding and Patterson function synthesis calculations were 
performed by the U. S. Naval Postgraduate School's CDC 1604 electronic 
computer. The diffraction intensity data was punched onto computer 
cards for input to a Fortran program developed by Prof. J. R. Clark, 
Dept. of Material Science and Chemistry, to execute the above calcu- 


lations. 


6. Experimental results. 

The hardness/ageing time curve for the specimens aged at 165°C 
1s shown in Fig. 1, page 9. The specimen aged at room temperature for 
65 days attained a Rockwell "F" hardness of 77.4. 

The lattice parameter of the as quenched specimen was 4.044A, 
which corresponds to 24 weight per cent copper in solution. The 
specimen aged at 165°C for one and two days retained this same lattice 


parameter. The lattice parameter of the specimens aged for five days 
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and longer was 4.0494 » which is the same as that of pure aluminum. 
The as quenched specimen retained the same lattice parameter on ageing 
at room temperature for 65 days. 

The Patterson functions for the artificially and naturally aged 
specimens are shown in Appendix I on both a small and large scale. 
The ordinate is P(w), where w is the interatomic vector length. The 
vector between any two atoms in the crystal always has its tail at 
the origin in the Patterson function. The abcissa then represents 
the magnitude of was it ranges from zero to "a", the lattice param— 
eter of aluminum, in increments of .Ola. The Patterson functions 
displayed in Appendix I are the one-dimensional projections of the 


electron density in the <\OO>,, direction. 


T. Discussion of results. 

The lattice parameter and hardness measurements were used to 
confirm the structures present at each ageing time. The lattice 
parameters were determined with only sufficient precision to detect 
any changes as a function of ageing time. The structures contained 
in each specimen were determined by comparing the experimental 
hardness curve with the hardness curves of Silcock, et ails It 
was not possible to obtain the complete sequence of intermediate 
precipitates on ageing at 165 C since only 24% copper went into 
solution. Initial ageing at 165°C produces G.P.\(2] and Q’ simul- 


taneously, while at peak hardness (five days) © predominates. The 


specimen fully aged at room temperature contains only G.P.\\] since 


SO Sanh ate et al., op. cit., p. 243. 
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the precipitation process does not progress beyond this point in 
naturally aged specimens. 

Lattice parameter measurements confirmed that the above struc- 
tures were present. A constant value of parameter is characteristic 
of the G.P. foneee It was observed that the parameter remained 
constant during natural ageing. The parameter also remained constant 
throvgh two days ageing at 165 C, indicating that G. P.(2] was present 
during this time. At five days ageing the parameter had relaxed to 
the pure aluminum value. Hence Q became the predominant intermediate 
at this time. 

When comparing the Patterson origin peaks, Appendix I, Figs. 4 
throvgh 11, it is observed that the as quenched and naturally aged 
specimens have peaks exactly at the origin while those of the arti- 
ficially aged specimens are shifted slightly. Also the origin peak 
for the naturally aged specimen is roughly twice the magnitude of 
all the others. It has been shown that the Patterson peak height 
is not a true measure of the product of the number of electrons, 
Z,22, contained in the two atoms separated by the vector that corres- 
ponds to the eek. * The peak height is directly proportional to 
Zi22 but inversely proportional to the sum of the cross-sections of 
the electron density peaks for a one-dimensional Patterson peak. 
Therefore a high Patterson peak is not necessarily due to two heavy 


atoms. The height also depends on the product of the electron 





lltpid. 
12 


A. Guinier, Solid State Physics, Vol. 9, p. 341. 


13Buerger, op. cit., p. 26. 
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densities of the atoms in the pair. 

The copper atom has a greater electron density than aluminun, 
Copper contains 29 electrons compared to 13 in aluminum and is 11% 
smaller in atomic size. Hence the presence of copper in an atom 
pair would substantially increase the Patterson peak height for that 
pair. In a specimen fully aged at room temperature a large fraction 
of the copper atoms are assembled in the G.P. zones, which are a few 
atomic layers thick and are on the average less than 100A apart. 14 
When looking perpendicular to the zones in the (\OO> , direction 
there would be a certain number of vectors of length "a" between 
atom pairs containing one or two copper atoms. There would be a 
smaller number of like vectors in the as quenched specimen, which 
is assumed to be a random substitutional solid solution of copper 
atoms in the aluminum matrix. The number of such vectors would also 
be less in the artificially aged specimens, which initially contained 
G.P.(2] and a) and later only 9. This can best be shown by consid- 
ering the accepted structure of 6, Fig. 2, page 13. The Q struc-— 
ture is tetragonal with parameters a = 4.04A and ¢ = 5.8A » and the 
axes are parallel to the C100 a) directions.!5 Only one of the 
three possible orientations of =) is shown in Fig. 2. When all three 
orientations are considered together, only two produce copper atom 
separations of 4.04A in the €\OO>, direction. Where G.P.(2] and 
ay occur together the copper atoms are shared between them and the 
total number available to share in atom pairs 4.04 A apart is less 

144. Guinier, Solid State Physics, p. 341. 


1593 Leock, et al., op cit., p. 242 
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than in the naturally aged specimen. This is also true when only 4 


is present. Hence one would expect G.P.(1t) zones to produce higher 
Patterson origin peaks than the other intermediate structures. 

As seen in Appendix I the left origin Patterson peaks occur at 
w values slightly less than "a" for all artificially aged specimens. 
This shift is real and is not due to an error in calculation or 
plotting of the Patterson function. The peak is at .98a for one days ’ 
ageing and at .99a for longer ageing times. This is best explained 
by considering the matrix distortion due to each of the intermediate 
precipitates. When a copper atom is substituted for an aluminum 
atom the matrix lattice contracts locally due to the smaller atomic 
size of copper. All three intermediate precipitates are coherent 
with the alwninum matrix. This coherency between matrix and preci- 
pitate lattices produces distortions in the matrix removed from the 
precipitate interface. The extent of the matrix distortion depends 
on the size and shape of the precipitate particle. Thus when G.P.[\\ 
zones are formed in their thin disc shape the lattice distortions 
are small and do not extend far from the interface. The average cell 
parameter is then that of the matrix and the Patterson peak remains 
at | w| = A. 

When G.P.(2)\ is present the lattice distortions are more extensive 
than those due to G.P.(\| , since the G.P(2) zones are thicker and have 
larger diameters. The distortions due to Q' are also extensive because 
of the 6 parameter of 5.8A. Thus when G.P.(2) and ey are present 
together after ageing one day the lattice is distorted by their 
combined effects and the Patterson peak occurs at \w | = .98a. AS 


i 
the G.P.L2\ zones dissolve and © grows with longer ageing times the 
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; 
lattice distortion decreases and is due to 9 alone at peak hardness 


of five days. 

In Appendix I, Figs. 12 through 19, are shown the Patterson 
functions on a large scale. In general, the number of resolved peaks 
increases as the ageing time increases. This is due to the increasing 
number of harmonics used in the synthesis of the Patterson function. 
The number of harmonics used in each case is determined by the point 
at which the coefficients of the broadened curve have decreased to a 
small value and those of the reference curve still have an appreciable 
value. The use of higher harmonics beyond this point results ina 
greater uncertainty in the coefficients of the corrected curve, 1® The 
intermediate peaks in the above Patterson functions are all quite small 
compared to the origin peaks. The resolution presented by the inter- 
mediate peaks is considerably greater than the 1.5 to 2A expected 
using x-radiation of 1.542A wave length. Thus the fine detail is 
false and is only a result of the harmonic analysis. No information 
concerning atomic positions can be gained from it. It is observed 
that the fine details of the Patterson functions for five days ageing 
and longer are essentially the same indicating that the.line profiles 


are very nearly the same. 


8. Conclusions and recommendations. 
The scattered intensity distribution in reciprocal space for 


the aluminum-copper alloy is shown in Fig. aes The scattering 





1©stokes, op. cit., p. 389. 


17,. Guinier, Solid State Physics, p. 339. 
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streaks around the nodes of the reciprocal lattice are due to the 
destruction of the strictly periodic electron density distribution 


in the crystal. This destruction is caused by the G. P. zones. 


002 








Fig. 3. Scattering in the reciprocal space of Al-Cu. 


The streaks pass precisely through the nodes of the matrix and along 
_the<\00> rows of the reciprocal lattice. The fact that the G.P. 
zones distort the matrix lattice causes the scattering streaks ton 
appear only on the high-angle side of the reciprocal lattice nodes. 
It was anticipated that the diffuse streaks about the nodes 
would contribute to the diffraction process in a distinctive way. It 
should be possible to differentiate G.P.(\) from G.P.[2] since the 
Mensth of the diffuse streaks acemceiinted to the size of the zone. 


This distinction could not be realized in the present investigation 


; \ 
as GeP.[2] zones were not obtained by themselves on ageing. However 
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it was possible to distinguish G.P.['} Peon Oere observing the 
height and shift of the Patterson origin peak. 

The amount contributed to the diffraction profile by the diffuse 
scattering streaks is uncertain. It is possible that the broadening 
of the profile due to instrumental factors, which were corrected for, 
is greater than that resulting from the diffuse streaks. However, 
Since the instrumental broadening remains constant with ageing the 
profile changes due to the diffuse scattering streaks could be detec-— 
table. The thermal vibrations of the atoms about their average position 
in the crystal also produce anisotropic diffuse diffraction effects 
about the nodes of the reciprocal lattice. The thermal diffuse 
scattering is related to the elastic properties of the crystal, which 
certainly change during ageing. Hence, the changes in thermal 
scattering should also SOntribUne to the line profile changes during 
ageing. The amount contributed and the best means of detecting it is 
a difficult problem to analyze quantitively. There has been no 
attempt in this investigation to do so. 

There are several ways in which the application of this method 
of analysis could be improved. The use of monochromatic x-rays wovld 
eliminate the profile broadening due to the spectral spread of the 
incident beam. In order to observe the effect of each intermediate 
precipitate on the line profile it is necessary to have them existing 
separately in the alloy. The transition in the aluninum-copper system 
from G.P.(\\ —»G.P.(2\ could best be observed by ageing a 4% copper 
alloy at 110 C, and from G.P.\2\—™ Q by ageing the same alloy at 190 C. 


The complete intermediate sequence could be achieved with definite 
separations in structures by ageing a 4.5% copper alloy at 130°C but 


would require ageing times of about 400 days. 
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APPENDIX I 


PATTERSON PROJECTIONS 


The Patterson functions for the naturally and artificially 
aged specimens are shown on a small scale in Figs. 4 through ll, 
and on a large scale in Figs. 12 throvgh 19. The ordinate is P(w), 
where w is the interatomic vector. The abcissa represents the 


magnitude of was it ranges from 0 to "a". 
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FIGURE 4. 
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Patterson Function, aged 65 days at room temperature 


FIGURE 6. 
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Patterson Function for specimen aged 1 day 


FIGURE 6. 
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Patterson Function for Specimen aged 2 days 


FIGURE 7. 











Patterson Functien for specimen aged 5 days 


FIGURE 8. 
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Patterson Functien for specimen aged 9 days 


FIGURE 9. 
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Patterson Functien for specimen aged 14 days 


FIGURE 10. 
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Patterson Function for specimen aged 30 days 


FIGURE ll. 
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Patterson Function for specimen aged 2 days 


FIGURE 15. 








Patterson Function for specimen aged 5 days 


FIGURE 16. 
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APPENDIX II 


FOURIER COEFFICIENTS FOR REFERENCE, BROADENED 
AND CORRECTED PROFILES 


The column headings for the tabulated Fourier coefficients 


are as follows: 


N, the harmonic number 
ALPHA, coefficient of the real component 
BETA, coefficient of the imaginary component 
INTENSITY for each harmonic is computed by 
ty = - Ours By ) 
where L is the number of equal intervals 
between the data points used to describe 
the curve. 
The corrected coefficients are normalized. The arrow in the 
"corrected" column indicates the number of harmonics used in the 


synthesis of the Patterson function. 
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30 DAYS, (200) 
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APPENDIX III 


EXPERIMENTAL DIFFRACTION LINE PROFILES 


The following instrumental settings were used in obtaining all 


of the diffraction line profiles. The abcissa is in degrees 29. 


Cu Kx radiation, 50 KV, 40 ma, Ni filter 

o 
Scan Speed % /minute, Chart Speed 30"/hour 
Time Constant 2 seconds, Scale 64 x .6 
Angular Aperture Tee Receiving Slit 0.003" 
Xenon Proportional Counter, 1725 volts 


Base 7 volts, Window 3 volts 
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APPENDIX IV 


THE FOLLOWING IS A FORTRAN 60 COMPUTER PROGRAM 
THAT WILL PERFORM THE STOKES UNFOLDING AND PATTERSON 
FUNCTION SYNTHESIS. 


PROGRAM UNFOLD 


THIS PROGRAM WILL COMPUTE 
le THE FOURIER COEFFICIENTS A(Ns91) AND B(Ns1) OF A CURVE 
2e THE FOURIER COEFFICIENTS A(Ns2) AND B(Ns2) OF A SECOND CURVE 
36. THE FOURIER COEFFICIENTS A(N) AND B(N) OF THE UNFOLDED CURVE 2/1 
4e THE INVERSE FOURIER TRANSFORM (SYNTHESIS) OF THE UNFOLDED CURVE 
5 THE CENTER OF GRAVITY OF THE UNFOLDED CURVE 
66 THE LOCATION OF THE MAXIMUM IN THE UNFOLDED CURVE 
7e A SERIES OF CURVES CAN BE SUCCESIVELY UNFOLDED 
THE RESULTS CAN BE PRESENTED IN BOTH TABULAR AND GRAPHICAL FORM 
GIVEN THE COEFFICIENTS IF(HsKsL)I**29 THE PROGRAM COMPUTES THE 
PATTERSON FUNCTIONs Y(UsVoW) (TeE« THE SUMMATION OVER THE COSINE 
‘TERMS OF A FOURIER TRANSFORM OF I(Fs(HsKsoL)I*¥*2 )e 


THE RESULTS ARE PRESENTED AS ONE OR TWO DIMENSIONAL SECTIONS 


ONE DIMENSIONAL SECTIONS CAN BE PRESENTED IN BOTH GRAPHICAL AND 
TABULAR FORM | 
TWO DIMENSIONAL SECTIONS ARE PRESENTED IN TABULAR FORM 
THE EXPERIMENTAL VALUESs IF(HsK9L)I*¥*2 ARE CORRECTED FOR LORENTZs 
POLARIZATION» MULTIPLICITY AND TEMPERATURE FACTORS 


FIRST DATA CARD CARRIES GENERAL TITLE - COLUMNS 2 = 80 


SECOND DATA CARD, . ; 
| COL 1-3: NA ='1 FOR THE FIRST CURVE OF A SET 
| . 2 FOR THE SECOND CURVE OF A SET 
QO IF FOLDING IS TO BE CALCULATED 
1 TO COMPUTE FOURIER COEFFICENTS ONLY 
A VALUE FOR NB MUST BE ENTERED OR STOP 15 
WILL RESULT. 
DIMENSIONS OF AAy BBy AND FIX DEPEND ON 
VALUE OF NA SPECIFIED 
COL 6-7 NC = +1 SYNTHESIS OF UNFOLDED DATA 
. = -1 SYNTHESIS OF NEW DATA 


. 


COL 4-5 NB 


me ce en re ee a eg ne re ee 





NOANDNNNDNNDNDANNDNDNNNNNDNNDNDNNDNDDNDNNNDNDNNNDDNDNONDNDNNNOONNNNNNANANGA 


COL 8-9 ND = 0 NO GRAPH OF PATTERSON FNe 
= 1 GRAPH PATTERSON FNe SEE BELOW. 


COL Gs 74 TITLE OF UNBROADENED CURVE 
THIRD DATA CARD = CONSTANTS FOR ANALYSIS OF UNBROADENED CURVE 
COL 2) =s%. iN = NO OF INTERVALS 
6 —- 10 NCG = -1 FOURIER SYNTHESIS =| GRAPH 


O NO SYNTHESIS. 
+1 FOURIER SYNTHESIS 
11 - 15 MODE = 1 IMPUT DATA PRINTED 
- = 0 IMPUT DATA NOT PRINTED 
16 - 20 J — FIRST OUTPUT HARMONIC 
| (NORMALLY = 0O) 
21 - 25 K - HARMONIC INCREMENT 
(NORMALLY = 1) 
26 - 30 L- FINAL HARMONIC 
| 31 - 35 K2- SUB HARMONIC —- IF Os PROGRAM 
Seis) Kk20 =) a : 
36 - 40 LC = 0 READS DATA CARD 
: MEMORY ° 
= 1 USES DATA AVAILABLE IN 
41 - 50 THO =INITIAL TWO THETA 
51 - 60 THF = FINAL TWO THETA 
61 - 73 XL = WAVELENGTH 


Won 


FOURTH DATA CARD —- UNBROADENED CURVE»N EQUALLY SPACED POINTS FORMAT 
(1l6F5e1l) » ON THE NECESSARY NUMBER OF CARDS« 


REPEAT SEQUENCE WITH A DUPLICTE GENERAL TITLE CARD » A SIXTH 
CARD SETTING NA = 2 AND CARRYING A TITLE FOR THE BROADENED 

CURVE » A SEVENTH CARD SETTING THE CONSTANTS FOR ANALYSIS OF THE 
BROADENED CURVE »AND AN EIGHTH CARD WITH THE N POINTS OF THE 
BROADENED CURVE. 


FOR A GRAPH OF THE PATTERSON FNs FOLLOW THE TWO SETS OF DATA CARDS 
(FOR THE REFERENCE AND BROADENED CURVES) WITH TWO MORE CARDS. 
THE FIRST ONE IS BLANK AND THE SECOND ONE SETS NA=NB=NC=05 ND=le 


e FURTHER CONVOLUTIONS CAN BE CARRIED OUT WITH 
ADDITIONAL SETS OF FOUR DATA CARDS » INCREASING NA BY ONE FOR 
EACH FOLD UP TO 10. 


ADDITIONAL CURVES CAN BE UNFOLDED BY ADDING ADDITIONAL SETS OF 
DATA CARDS SETTING NA =1 AND 2 IN THE SECOND AND SIXTH CARDS OF 
BACH SET 


FOR A FOURIER ANALYSIS SET NB = 1 IN COL 5 IN THE 
SECOND CARD OF A SET OF DATA CARDS 


ee Sant 
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NA ANAND ANDNANANNAAADAAANDADNDA 


THE PROGRAM IS TERMINATED BY TWO BLANK DATA CARDS 


STATEMENT NUMBERS 
m2 3 4 S 69°F) 8 9 10°91 12°81 ea cee 
18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 
35 36 37 38 39 40 41 42 43 44 45 46 
51 52 53 | 
60 61 


Vx 
i 
/ 


: 103 104 105 106 107 108 109 
io 11) 112 113 114 115 116 117 118 119 
ie 121 122 123 124 125 126 127 128 129 
130 131 132 133 134 135 136 137 138 139 
140 141 142 143 144 145 146 147 148 149 
eels] 1525153 154 155 156 1s7elee 159 
160° 161 162 163 164 165 166 167 
501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 
900 901 902 904 906 907 908 909 
orea48 49 70 170 71 72 73 


1 FORMAT (1H1////) 
2 FORMAT (80HO 
i ) 
24 FORMAT (159 2F10e4s Fl005 9 2F1004 sF1l0e5 s2F1l004s F1l005) 
39 FORMAT(I155 2F15e4s Fl5e4//) 
40 FORMAT(1HO13X1HN1LOXSHALPHALIX4HBETA6X QSHINTENSITY) 
900 FORMAT ( 139 312s 8A8) 
901 FORMAT( 81559 3F10e3) - 
902 FORMAT (16F5el) . 
904 FORMAT(3HON=14,5 #£=°7H J=155 7H K=15s 7H L=15+¢ 
> 1 8H K2=1597X6HNCG = 13///12H DATA POINTS//) 
906 FORMAT (1HO 16X SHREFERENCE 20X 11HBROADENED 21X QHCORRECTED //) 
907 FORMAT(1HO 4X 1HN 2X SHALPHA 7X 4GHBETA 3X SHINTENSITY 4X SHALPHA 
-  1,6X 4HBETA 2X 9HINTENSITY 4X SHALPHA 6X 4HBETA 2X 9HINTENSITY//) 
908 FORMAT(1X) 
909 FORMAT(1Xs 16F5el) 


DIMENSION Y¥(200)s AA(20094)s BB(20094)s FIX(20094)s NAME(8B) » 
1 FR(20094)5 FI(20094)s FIT(20094)s 
2 21919900) . 


ian J 
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ST 
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al 


31 
30 


23 
2 


oa 
Bo 13 


1 
514 
512 
13 


: H , 38 


DIMENSION WW(900)»)  ITITLE(12) 
Bon 141i 


IVitEEt 1) (= 8h 
ITITLE(1) = 8H LARSEN 
TT TTLE C2 eA emesH 0276 
ITITLECT)° = BHPATTERSO 
ITITLE(8) = 8HN FUNCT 
PTITLE( 9) = SHION: JFOR 


ITITLE(10) = 8H AL-CU e 
LABEL = 4H 65D | 
TPL = 6028318530718 

R=100/(57e2957795131%200) 

READ 2 

READ 900» NAs NBs NC» ND» NAME 

MeNGraboes. Ib. 47 

IF(NA) 155 159 9 . 

IF (ND) 48948949 . 


CALL DRAW(1019Zs9WWs090 sLABEL s ITITLE 909100090 9090909991091 LAST) 


GO TO 3 
STOP 48 
READ 901» Ns» NCGs MODE» J» Ks Ls K2e9 LCs THOs THF» XL 
PatkK2)} ~30551%530 S 


N 


ee= 1 
FN = N 
TNR 20/FN 


THF = THF#R 
THO = THO * R 
IF LC EQUALS ZERO» READ NEW DATAs OTHERWISE OPERATE ON OLD 
ee (LG W535 2633932 

READ 9029 (Y(I)s T= lsN ) 


ne 
. 


PRINT 9005 NAs NBs NCs NDs NAME 

IF (MODE) 125135312 F 

PRINT 9049 No Js Ks Lo K2s NCG 

PRINT 909, (Y{E)> TAa aie) 2 
IF (NA —- 1) 5155 5139 514 
STOP 515 

AREA = 0.0 

DO 511 I = lo»N | 

AREAL =(Y(1)*#(SINF( THF) = SINF(THO)))/7 (XL * FN) 

AREA = AREA + AREAI 

DO 512 I = 1oN 

Y (Ll )v=a¥ ite 7 AREA 

XK2 = K2 

XJ = J 

ARG = TPI/FN/XK2 
ARG1 = ARG*#¥XJ 


t 
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CS = COSF(ARGI) 
TCS = 2e*CS 


SS = SINF(ARG1) 
Kak 

ARG2 = ARG*XK 
CB = COSF(ARG2) 
SB = SINF(ARG2) 
M1l=J 


GO TO (50155025503 3504 35055506450/550655093510)45 NA ya 
501 NAN=NA 
- GO TO 14 S 
502 NCW=NA | 
NAN=NA=-1 
GO TO 14 
503 NCW=NA-1 
NAN=NA=2 
| GO TO 14 
504 NCW=NA=2 
NAN=NA=3 
GO TO 14 
505 NCW=NA=3 
~~ NAN=NA=4 
GO TO 14 
506 NCW=NA-4 
NAN=NA=5 
GO TO 14 
507 NCW=NA-5 
) NAN=NA=6 
GO TO 14 
508 NCW=NA—6 
NAN=NA-7 
GO TO 14 
509 NCW=NA=-7 
NAN=NA-8 
GO TO 14 
510 NCW=NA—8 
3 NAN=NA=9 7 
BEGINNING OF LOOP ON Ml 
14 T3 De 
Oe 
N 
18 ‘Ti Tess 12 =.13 #2714 
I = [=] 
Pesto 16s | Sen 7 


—i 
tot ott 





20) 
al 


37 


43 


38 


28 


Z2 


45 
46 


35 


Ne RS LR ET A A RT A, SL, KR AR, OS EN fe La a hg Se nr ae ge Sa =< ome see is ~~: fo on, = 
: ' 


AA(MsNA) (T1\- CS*T2)*TNR 
BB (M»sNA) (SS* T2)*TNR 
FIX(MsNA)=(AA(M9NA)*¥AA(MsNA) + BB(MyNA) * BB(MaNA))/TNR 
TEST FOR END OF Ml SERIES 

IF (M1l-L) 21526520 

STOP 20 

Ml=M1+K | | 

TCS = CB¥ES = SB¥SS 

CB*SS + SB*CS ‘ 


IF (NO) 239425 27 
STOP 23 

IF (NB) 2793927 
PRINT 1 

PRINT 2 

L2=L+l1 

J2=Jtl1 

IF(NB) 37938937 
PRINT 40 

PRINT 908 

DO 43 M = J2s L2 
M2=M-1 
PRINT 399M2sAA(MsNA) sBB(MoNA) 9FIX(MsNA) 
GO TO 3 

PRINT 906 

PRINT 907 

IF (NA=2) 22522928 


Ak 
“~ 


AA (MsyNAN) =AA(MyNO) 

BB (MsNAN) =B8B(MsNO) 

FIX(MyNAN) = (AA(MsNAN)*AA(MsNAN) + BB(Ms NAN) *BB (Ms NAN) )/TNR 

DO 25 M=J29L2 

M2=M-1 

J2=J +1 

FR(M»NO) = (AA(MyNCW)* AA(MsNAN)+ BB(MsNCW)*BB (My NAN) ) /FIX(MyNAN) 
1 *TNR : 

FI(MyNO) = (BB(MyNCW)* AA(MsNAN)— AA(MsyNCW)* BB (Ms NAN) )/FIX(MyNAN) 
1 ¥TNR : — 


IF(M -J2)465 455 46 
FA = FR(J29NO0) 


FR(MsNO) = FR(MsNO)/FA 
FI(MsNO) = FI(MsNO)/FA 
FIT(MsNO) =(FR(MsNO)*FR(MsNO) + FI(MsNO)* FI(MsNO))/TNR 


PRINT 24 sM2sAA(MsNAN) sBB(MsNAN) s9FIX(MsNAN) 9 AA(MsNCW) 9 BB( Ms NCW 


1) » FIX(MsNCW) »sFR (MseNO) » FI(MsNO) sFIT(MsNO) 


TF (XMODF(M 945))29935929 
PRINT 1 


‘ 


\ 
60 | 





te ee 


Camara are a@raG) Gionts @y OO) OY CY GO) OVAT) OO. AAD OAROANANNNNOADANG 


29 
36 


Za 


BRENT 2 
PRINT 906 | 
PRINT 907 
BO, 10425 

IF (XMODF(M 
PRINT 908 
Ge! TO 25 
CONTINUE 
i tNCG) 


¢9))2595 36925 


1225 44s, 122 


PROGRAM FOUSYN 


1ST DATA CARD 
_ COL 1-80 TITLE 


2ND DATA CARD 
pCOL Jie— 4 DEBYE TEMP. 


B 


B=0.0 IF TEMPe FACTOR IS NOT COMPUTED 


Come it LAMBDA 


- 11 = 14 NUMBER OF 
-COL 15 - 20 NUMBER OF 
21 - 25 NUMBER OF 
26 - 30 NUMBER OF 


rEAKS USED 


INTERVAES/ CEE IN X DIREGihReNS 
INTERVALS/CELL IN Y DIRECTION 
INTERVALS/CELL IN Z DIRECTION 


ADJUST DIMENSION STATEMENT. FOR FRACTIONS UVW OF XYZ USED 


3RD DATA CARD 
COL 


eee MAX VALUE OF TU 


4 —- 6 MAX VALUc OF V 


7 - 9 MAX VALUE OD W 
10 - 12 MAX VALUE 
13 - 15 MAX VALUE 
16 - 18 MAX VALUE 

| 19 - 21 MAX VALUE 
21 - 23 MAX VALUE 
24 - 26 MAX VALUE 
WANTED | 


ADJUST DIMENSION STATEMENT 


4TH DATA CARD 


COL 1-3 H 
rt am OK 
7T-9OL 
10 - 
“15 = 18 TWO THETA 


ieee OeMUE | 1 Pie C llny 


% 


eZ 


U FOR PRINTOUT 
V FOR PRINGOUR 


FOR MAX VALUES OF HKL USED 


14 RELATIVE INTENSITY 


6l 


W FOR PRINTOUT ACCORDING TO 


(WITH DECIMAL POINT) 


\ 


« 
A eS A SS SS en eS i ES rn, Se 


SECT LONG 


— ey ee ee oe 7 — — 





122 CONTINUE 


51 FORMAT (3(3135 F5eOs F4e09 129 F5e0)) 
52 FORMAT(F4.e2sFle5s 145 315 ) 
53 FORMAT (6139 3F4e¢2s5 2F10e4) 
4 FORMAT (1HO»s2X»s 1HNs 6Xs 22H( Ue Vii W )9Xs°9H Y(UsVow) 
l Le) 
5 FORMAT (155s 6X9 F604 9 1Xo F6e4s 1X9 Feds 6Xo F1204) 
6 FORMAT( 50H 
. 7 FORMAT (1X) 
8 FORMAT (1H1//// ) 
O FORMAT(21HOCENTER OF GRAVITY = F602) 
1 FORMAT (55HOMAX HARMONIC USED FOR SYNTHESIS OF CORRECTED CURVE = 
eo 3/7) : 


167 DIMENSION F(1l919500)s TTH(1919500)s MU(1919500)» C3(100)> 
i 23 (VCO IC] 5:19 20 Ones Ciel 1 core 2 (1) 511 aoc 


| TF(NC) 129) 166% 128 

166 STOP 166 

129 READ 6 
READ 529BsXLsNP»NUsNV s NW 
READ 539 Jl» J2» J3»9 IT1MAX» I2MAX»s U3MAX» Ul» Vly Wls TTHOs TTHF 
GO TO 136 . 

128 NN1 = 0 


ass 
a 


0 
0 
0 


&THE leas 
136 IT1M = T1MAX + 1 
IT2M=I 2MAX+1 








13M I3MAX +1 
TTHD TTHO - TTHF 
Pow lls [32 lon 
1S3 = (13 = 1 e1oo 
DO 115 12 3 ls 12M 
1S2 = (12-1)*10 + 153 
PO 115 11=1,11M 
1S=1S2+11 
F(I]S)=0e0 
FIC(IS) = 020 
é TPM is) =0.0 - 
peo) Z7(1S)=0.0 
hem NG) F275, 155 126 
mes DO 137 M = 1» L2 
IF (M=1) 159391595160 


159 I4M =0 

160 IF(I4M) 16151619158 

161 F(ls ls M) = FR(Ms NO) : 
Fl( ls ls M) = FI(Ms NO) ‘ 
MCE RI Ms NOI I 17050 1705 13g 

170° I4M = M 
SO 10 137 

158 F(ls ls» M) = FR(Ms NO) 


Fl(ls ls M) = FI(Ms NO) 
137 CONTINUE - 
FOR TRY = -1e05 Y(UsVseW) IS SYNTHESIZED USING THE FIRST M° 
HARMONICS THAT KEEP Y(UsVoeW) POSITIVEe~ 
Mowe ik = -1.20 . 
13M = 14M -] 
ISM = I4M -2 
SO 10 125 ' 
FOR TRY = 120% Y(UsVeW) IS SYNTHESIZED USING ALL L HARMONICS. 


ee, TRY 10 


Mogmir (XMODFINP, 3))111s 112» 121 

112 NP = NP/3 
GorTO 113 

111 NP = (NP/3) +1 

113 DO 114 NR = l» NP | 
READ 51+ NNl»s NN2»s NN3s FLIS)» TIH( ISI» MULISIs F1C15) 
1S=NN3¥100+NN2¥104NN141 
THA TS) =1 Ht 1S) eR 
AAA=SINF(TTH(1IS)) 
AAA2=AAAK# 2 | 
TE=1 00 /EXPF (B¥AAA2/XL¥¥2) 

- - 2Z=COSF(2.0*TTH(IS)) 





A eR, A ne NG fe 


POL=(1-0+Z2Z**2)/AAA2Z*COSF(TTH( ES) } 


AMUCTS) = MU(ITS)} 
PLCIS) =F L01S) /( Tee Ome Mut IS)) 
m4 F(1S) = F(iS)/(TE * POE a aMutiS)) 


125 PRINT 8 
IF (NC) 1315 155 130 
131 PRINT 6 
GO TO 132 
130 PRINT 2 
PRINT 61s I5M 
132 PRINT 4 
TP1=6.28318530718 
MARK=0 
IF(NC) 157, 138, 139 
138 STOP 138 


mo JU = 1 

a= i 

GO TO 140 | 
[e700 104 JU= 1, J1] 

De oe JY Sy : 
140 DO 104 UW = ly» y33 


eeelc (NC) 153551545 155 
154 STOP 154 
no |) = 


153 AJU 


V = AJV/ANV oa 
156 AJW = JW - 1 oe 
ANW = J3 
wo= AJW/ANW 
Pe ONG jess) 142, 142 
142 STOP 142 


f 


143 11 = 1 
P=. 
GO TO 144 


141 DO 103 I1=1,11M 
DO 103 12=1,12M 
144 DO 103 13=1,13M 


NN1 = 11 - 1 
NN2 = I2= 1 
NN3 = 13 - 1 
FNL = NN1 
FN2 = NN2 
FN3 = NN3 





fe 


TS RI en ee 





IF (NC) 14559 1465 147 


fey Sl1{(11) = 1.0 
S62 12) 41.0 
GO TO 148 

145 S1l(I1) = SINF(U*FN1¥*TPI ) 
S2(12) = SINF(V¥FN2*TPI) 

148 S3(13) = SINF(W*FN3*TPI) 
TS = Sl(I1) * $2012) * S313) * FEC I1912913) Fa 


be ING), 149i a O's 15 1 


eo Clt( ll) = GCOSF(UFFNI1*TPI) 
G2) 2) COSF(V#¥FN2*TPI ) 

P52 €3(13)=COSF (WeFN3*TPI) 
B= Gl Gin G22) eC3 ( 13) | 2 ol3) 
me 1S Wiel ys 6116 D7 

116 TSM=0.0 

mea etSM = TSM + T + TS 

103 Z(JUsJVsIW)=TSM 

107 MARK = MARK +1 
JJW = JW - 1 | 
PRINT 59 JSJIWs UsVoWs Z(JUs JIVsJIW) 
WW( JW) = JJW 
IF (XMODF (MARK»s 35)) 1095108:109 

TOS [F (NC) 134s 15,5 133 

134 PRINT 8 
PRINT 6 
GO TO 135 

M35 PRINT 8 
PRINT 2 

me> PRINT & 
SO 10 165 . 

109 IF (XMODF (MARK » 5)) 16591109165 

110 PRINT 7 

toon CONTINUE 

104 CONTINUE 

| IF (TRY) 72944944 

44 (G60 TO 3 

72 ‘IF (ND) 16391639 73. 

73 CALL DRAW(1019Z sWWs0s0sLABELsITITLEs09100905 0, 050+ 6s 1 Onn LAST) 
GO TO 163 : 
END 
END 


f ae 








—— 




















